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ABSTRACT 

The recent detection of a 3-hr X-ray flare by the Chandra Observatory has 
raised the possibility of enhanced emission over a broad range of wavelengths from 
Sgr A*, the suspected 2.6 x 10 6 M black hole at the Galactic Center, during a 
flaring event. We have, therefore, reconstructed 3-hr data sets from 2/im speckle 
and adaptive optics images (9 core = 50 - 100 mas) obtained with the W. M. 
Keck 10-m telescopes between 1995 and 2001. In 25 separate observations, no 
evidence of any significant excess emission associated with Sgr A* was detected. 
The lowest of our detection limits gives an observed limit for the quiescent state 
of Sgr A* of 0.09±0.005 mJy, or, equivalently, a dereddened value of 2.0±0.1 mJy, 
which is a factor of 2 lower than the best previously published quiescent value. 
Under the assumption that there are random 3-hr flares producing both enhanced 
X-ray and near-infrared emission, our highest limit constrains the variable state 
of Sgr A* to <0.8 mJy (observed) or 19 mJy (dereddened). These results suggest 
that the model favored by Markoff et al. (2001), in which the flare is produced 
through local heating of relativistic particles surrounding Sgr A* (e.g., a sudden 
magnetic reconnection event), is unlikely, because it predicts peak 2//m emission 
of ~300 mJy, well above our detection limit. 

Subject headings: Galaxy: center — infrared: galaxies — X-rays: galaxies - 
accretion, accretion disks — galaxies: jets — black hole physics 

1. Introduction 

The variability of Sagittarius (Sgr) A* at X-ray wavelengths (Baganoff et al. 2001a) has 
bolstered the case for associating this source with the suspected 2.6 x 10 6 M black hole at 
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the center of our Galaxy (Eckart & Genzel 1997; Genzel et al. 1997, 2000; Ghez et al. 1998, 
2000). In two Chandra observations separated by almost a year and having a total of 76 ksec 
of exposure time, Sgr A* was detected at X-ray wavelengths for the first time and was also 
seen to flare in intensity over a time scale of 3 hours (Baganoff et al. 2001a,b). While the 
flare's short duration implied a small region of origin, <400 R s (where R s is the Schwarzschild 
radius = 2GM./c 2 ), its large amplitude, a factor of 50, has raised the possibility of detecting 
corresponding intensity enhancements at wavelengths outside the X-ray regime. 

Existing models for Sgr A*'s flared state make very disparate predictions for the emission 
at wavelengths between the X-ray and radio regimes (Markoff et al. 2001; Liu & Melia 2002). 
The wide differences between these models are a result of assuming different geometries (disk 
vs. jet) and emission mechanisms for the flaring process (e.g., enhanced accretion rates vs. 
magnetic reconnection). In some models, the predicted emission in the flared state, at 
infrared (IR) wavelengths, dramatically exceeds that of existing detection limits (Genzel & 
Eckart 1999; Stolovy et al. 1999; Morris et al. 2001). For example, the preferred model of 
Markoff et al. (2001) predicts an observed 2^m flux density of ~13 mJy, or a dereddened 
flux density of ~300 mJy, during the flared state. Unlike the situation at radio wavelengths, 
where long-term monitoring campaigns have been used to constrain the flared state of Sgr A* 
(Zhao, Bower, & Goss 2001), the limited time coverage and spatial resolution of published 
IR experiments prevent meaningful constraints on the flared state's IR emission from being 
inferred 4 and, thus, the reported limits are assumed to be associated with Sgr A*'s quiescent 
state. 

The W. M. Keck Observatory dynamical study of stars in the central stellar cluster 
(Ghez et al. 1998, 2000; Gezari et al. 2002) provides a rich source of high angular resolution 
2/im data between 1995 and 2001. In this paper, we present 2//m flux density limits from 
maps that were each composed of data from a single night. The elapsed time of 3-4 hours in 
each map is approximately the same as the time scale of the observed X-ray flare, making 
this data set ideally suited for possibly detecting a flare of this type. Given the quantity of 
these observations, our non-detections establish a robust upper limit for the flared state's 
2|Um emission intensity. 



4 While several papers have reported the possible detection of a variable near-infrared source coincident 
with Sgr A* (Hcrbst, Bcckwith, & Shurc 1993; Close, McCarthy, & Melia 1995; Genzel et al. 1997), subsequent 
high resolution observations have identified this emission to be from high proper motion sources (Eckart et 
al. 1995; Eckart & Genzel 1997; Ghez et al. 1998, 2002). 
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2. Observations 

High resolution, near-infrared observations of the Galactic Center were conducted from 
1995 June to 2001 July using both speckle and adaptive optics (AO) imaging techniques 
on the Keck 10-m telescopes. The speckle observations were obtained in the K-band (A c 
= 2.2/iin, AA=0.4yum) using the Keck I facility near-infrared camera (NIRC; Matthews & 
Soifer 1994) with external re-imaging optics. This resulted in a pixel scale of 0"0203 and a 
field of view (FOV) of 5 / . / 12x5'. / 12 (Matthews et al. 1996). During each night of observations 
several thousand short exposures (t exp — 0.137 sec) were taken in sets of ~200. A more 
limited set of data was collected using two different science cameras behind the Keck II AO 
system (Wizinowich 2000). The first AO data set was collected in the K'-band (A = 2.1/xm, 
AA=0.35yum) in early 1999 with the near-infrared engineering camera (KCAM; Wizinowich 
et al. 2000), which had a pixel scale of 0"0175, a FOV of 4"4x4"4. Each image had an 
exposure time of 5 sec. The slit-viewing camera of NIRSPEC (SCAM; McLean et al. 1998) 
provided a second set of AO images for this study. These images, like the speckle images, 
were made in the K-band and have a pixel scale of 0"0170 and a FOV of 4'. / 4x4 / .'4, and 
exposure time of 10 sec. USNO 0600-28579500 served as the natural guide star for all of 
these AO observations. Since this guide star is both faint (R = 13.2) and distant from the 
target (r~30"), the AO performance was non-optimal. Table 1 provides a summary of all 
observations. 



3. Data Analysis & Results 

Three basic steps constitute the data analysis process in this program. First, high 
angular resolution maps are generated from the individual short exposure frames (§3.1). 
Second, all point sources in the FOV are identified and a direct detection of Sgr A* is ruled 
out (§3.2). Third, limits for Sgr A* are derived from the residual maps, in which all identified 
point sources have been removed (§3.3). 

3.1. Construction of Images 

Image processing proceeds similarly to that carried out for the dynamical experiment 
with one exception. Rather than combining all the data from the duration of an observing 
run, typically 2-3 nights, we synthesize the data over each night to produce 27 maps, each 
of which is limited to an elapsed time of 3-4 hours. Since the details of this method are 
described elsewhere (Ghez et al. 1998), only a brief summary is provided here. Standard 
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image reduction techniques are applied to all the individual speckle and AO frames. For 
the speckle data, a two stage shift-and-add (SAA; Christou 1991; Ghez et al. 1998) analysis 
then produces the final high resolution maps. In the first stage, the 200 frames in each set 
are combined to form an intermediate SAA image. Then, these multiple intermediate SAA 
images (from throughout the night) are combined to form one final SAA map. This allows 
each intermediate image to be examined for seeing quality. In combining the intermediate 
images, a seeing cut is established so as to exclude those images with the worst seeing from 
the final map. For the AO data, this cut is also carried out on the individual AO images 
before they are registered and averaged together. Figure 1 displays representative final 
nightly speckle and AO maps. 

3.2. Point Source Identification & Search for Sgr A*'s near-infrared emission 

In all maps, stars are identified using StarFinder, an IDL package developed for astrom- 
etry and photometry in crowded stellar fields (Diolaiti et al. 2000). This package iteratively 
generates estimates of the point spread function (PSF) from a few selected bright stars and 
then identifies point sources over the entire FOV through cross-correlation of the map with 
the PSF model. For the PSF extraction, we found that the most reliable PSF models are 
obtained with a support size of ~2", which represents a compromise between needing to 
accommodate the large PSF halos and yet having a limited FOV. This choice limits our 
analysis to images with PSF halo sizes of 0V4 or less, as the PSFs of the remaining two 
images are poorly characterized by this process. The PSF model is based on four of the five 
brightest stars in the FOV (IRS 16NE, 16C, 16NW, and 29N; see Figure 1); IRS 16SW is 
avoided as a PSF model star, since it is surrounded by relatively bright stars in its immediate 
vicinity. For point source identification, a correlation coefficient greater than 0.8 between 
the PSF model and the actual stellar image is required to avoid spurious detections. This 
process results in the identification of ~100 point sources in each map. 

The speckle and AO images have significantly different PSFs. Nonetheless, both PSFs 
are composed of a compact core on top of a broader halo. Table 1 provides the characteriza- 
tion of the PSF in each map based on the radial profile of the PSF model. While the speckle 
images have PSF core FWHM that are nearly diffraction limited (~0'.'05) and ~40% smaller 
than that of the AO images (~0'.'08), the AO PSF core contains ~30% of the total energy, 
~12 times more than the typical speckle PSF. 

At this stage in the analysis, it is possible to look for direct detections of Sgr A*. We 
use proper motion acceleration vectors (Ghez et al. 2000) to pinpoint the central black hole's 
position relative to the nominal radio position of Sgr A* (Menten et al. 1997) to within 0'.'04 
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(lcr). Within 0'.'08 of this location, 4 sources (13.9 < K < 16.5) are identified, all of which 
were previously detected in "monthly" maps made from all data in a single observing run 
(Ghez et al. 1998, 2002) and, furthermore, have significant proper motions. This high stellar 
density emphasizes the need for improved accuracy in Sgr A*'s position in the IR reference 
frame in order to measure or constrain its emission. With no stationary source identified in 
this region, we conclude that Sgr A* has not been detected. 

3.3. Flux Density Limits for Sgr A* 

In order to determine an accurate detection limit at the position of Sgr A*, it is necessary 
to remove the contaminating seeing halos from nearby sources. A 'stars only' map is created 
using the PSF model and list of stars generated by StarFinder. This is then subtracted 
from the original map, producing a residual map. With the residual map, a 3a point source 
detection limit for Sgr A* is established based on three times the RMS of 25 aperture 
photometry values, which are calculated using ~60 mas radius apertures and sky annuli 
extending from ~60 mas to ~90 mas. The 5x5 grid of apertures in the residual map 
corresponds to an area of ~0 / .'6x0 / .'6, more than two orders of magnitude larger than the 
uncertainty in the location of Sgr A*. 

Zero points are obtained through carrying out the same aperture photometry in the 
original maps (prior to the 'stars only' subtraction) on all known non- variable sources brighter 
than K=10.5,using the flux densities reported in Blum, Sellgren, & Depoy (1996), and that 
occur in more than 30% of the frames for each night. The photometric calibration sources 
used are IRS 16NW, 16C, 16CC, and 16NE, when the FOV allows its inclusion; IRS 29N is 
omitted as it is found to be marginally variable at the 2a level. Typical photometric zero 
point lcr uncertainties of ~0.04 mag result from this procedure. 

Table 1 and Figure 2 contain the resulting 3a point source detection limits for Sgr A*. 
The lowest of these upper limits gives an observed limit for the quiescent state of Sgr A* 
of 0.09±0.005 mJy, or, equivalently, a dereddened value of 2.0±0.1 mJy, while the highest 
limit constrains our analysis of the variable state of Sgr A* to <0.8 mJy (observed) or 19 
mJy (dereddened). Although these upper limits vary significantly from epoch to epoch (due 
to variable observing conditions and/or a variation in the number of contaminating sources 
detected and removed) the lowest of them is lower than the best previously reported limits 
at 2(Um for Sgr A*'s quiescent state (dereddened 4 mJy; Genzel & Eckart 1999) by a factor 
of 2. 



- 6- 



4. Discussion 

Using the coverage of the X-ray and near-infrared experiments and assuming random 3- 
hr flaring events that produce both enhanced X-ray and near-infrared emission, we consider 
the likelihood that a near-infrared flare in excess of our weakest limit occurred over the 
course of the IR experiment. The probability of seeing such a flare is given by the following 
binomial distribution: 

n\ 

P(u detections in n trials) = —. rrP U( l n u (1) 

v\\n — v)\ 

where p is the probability of detecting a flare in one trial and q is the probability of no 
detection (1-p). Here we will refer to p as the duty cycle and note that in our context, 
this quantity describes what fraction of a set of 3-hr images should contain a flare. Figure 
3 shows the probability distributions as a function of underlying duty cycle for both the 
X-ray experiment, in which v—\ and n=7.6 (obtained by dividing the total time, 76 ksec, by 
the length of the observed X-ray flare, 10 ksec) and the near-infrared experiment, in which 
u=0 and n=25. The X-ray and near-infrared probability distributions have only a small 
overlap, which suggests that if the flaring activity at X-ray and near-infrared wavelengths 
are coupled, a flare most likely occurred during the near-infrared experiment. The joint 
probability distribution quantifies this and suggests there is, at most, a probability of 9% of 
a random 3-hr near-infrared flare in excess of 19 mJy (dereddened) . We therefore assume 
that at the 2a confidence level, a flare occurred during our experiment and use our limits to 
constrain the variability models. 

While only a limited amount of modeling of the recent 3-hr X-ray flare detected at Sgr 
A* has been carried out, existing models predict 2^/m dereddened emission as high as 300 
mJy in the model preferred by Markoff et al. (2001) but as small as 0.4 mJy in Liu & Melia 
(2002). The former model explains the elevated X-ray emission, produced by synchrotron 
self-Compton, by an enhanced temperature for the relativistic electron population, as might 
arise in a magnetic reconnection event. On the other hand, Liu & Melia present a flare model 
in which the flare arises due to an enhanced accretion rate and bremsstrahlung emission is 
dominant at both near-infrared and X-ray wavelengths. The lack of a near-infrared detection 
of Sgr A* makes the Markoff et al. (2001) model and any other mechanism that produces 
flared 2/im emission in excess of 19 mJy (dereddened) unlikely. 

5. Conclusions 

This paper summarizes a search for a near-infrared counterpart to Sgr A* in the flared 
state. From the length of our observations, this search was sensitive to variability on time 
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scales of 3 hours. No such counterpart was detected. However, by identifying and removing 
all the stars in the crowded inner ~0'/6x0 / . / 6 of the Galactic Center, an upper limit for the 
emission from Sgr A* has been inferred for each observation epoch. These limits constrain 
the quiescent emission from Sgr A* to <0.09 mJy (2.0 mJy, dereddened) and the variable 
component to <0.8 mJy (19 mJy, dereddened) at the 2a confidence level. More X-ray data 
will improve the estimated flaring duty cycle and is likely to increase the confidence level of 
our limits for the near-infrared component of the flared emission. 

Valuable input on the manuscript was contributed by Gaspard Duchene, Caer McCabe, 
and Andreas Eckart. Data presented herein were obtained at the W.M. Keck Observatory, 
which is operated as a scientific partnership among the California Institute of Technology, 
the University of California and the National Aeronautics and Space Administration. The 
Observatory was made possible by the generous financial support of the W.M. Keck Foun- 
dation. This work was supported by the NSF both through grant No. 9988397 and the 
Science and Technology Center for Adaptive Optics, managed by the University of Califor- 
nia at Santa Cruz under cooperative agreement No. AST - 9876783. AMG also thanks the 
Packard Foundation for its support. 
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Table 1. Near- Infrared Limits on Sagittarius A 
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a Percentage of total energy contained in the PSF core 

b Column 7 lists observed limits and column 8 lists dereddened values using an A v =30 with Aj,/A„=0.112 and Ay /A t) =0.117 
(Melia & Falcke 2001; Rieke & Lebofsky 1985) 

c Limits for the two images with halo FWHM >0 / .'5 are not included in the analysis (see §3.2). 
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Fig. I. — Speckle (left) and adaptive optics (right) images from May 1999 of the Galactic 
Center. The small box indicates a l"xl" region centered on Sgr A*, whose approximate 
position is marked with a cross. Both images are displayed with a histogram equalization 
stretch to show the fainter stars in the field and are oriented such that North is up and East 
is to the left. 

Fig. 2. — The 3a limiting flux density calculated for Sgr A* for each epoch of observation 
(corrected for reddening by a factor of ~22.) The lowest limit of 2.0 mJy, a factor of two 
lower than previously published limits, and the highest limit of 19 mJy, constrains Sgr A*'s 
quiescent and variable states, respectively. 

Fig. 3. — The probability (P n ) of seeing the number of detected flares (v) of both this study 
(0 detections) and the Chandra experiment (1 detection; Baganoff et al. 2001a) for various 
duty cycles. The low joint probability suggests that if the flaring activity at near-infrared 
and X-ray wavelengths are coupled, a flare occurred during the near-infrared experiment and 
is restricted to be lower than our detection limit. 
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